UW/PT 00-01 

Sphaleron rate in the symmetric electroweak phase 



Guy D. Moore 

Department of Physics, University of Washington, Seattle WA 98195-1560 USA 

(January 2000) 



Abstract 



Recently Bodeker has presented an effective infrared theory for the dynamics 
of Yang-Mills theory, suitable for studying the rate of baryon number violation 
in the early universe. We extend his theory to include Higgs fields, and study 
how much the Higgs affects the baryon number violation rate in the symmetric 
phase, at the phase coexistence temperature of a first order electroweak phase 
transition. The rate is about 20% smaller than in pure Yang-Mills theory. 
We also analyze the sphaleron rate in the analytic crossover regime. Our 
treatment relies on the ergodicity conjecture for 3-D scalar c/)^ theory. 



Baryon number is violated in the standard model , and while the rate of the violation 
is negligibly small in vacuum, it can be significant at higher temperatures 0,^- 

Recently there has been substantial interest in how fast baryon number is violated in 
the standard model at high temperatures. This is both because the rate is important for 
baryogenesis, and because it is an interesting problem in the more general framework of 
understanding dynamics of hot gauge theories. 

The last 3 years have seen major progress in this problem. First, analytic work has 
clarified what physics is relevant. While it has long been known that infrared or "soft" gauge 
fields with momenta p ~ g^T are responsible for baryon number violation 0, the recent 
work of Arnold, Son, and Yaffe has shown that "hard" modes with momenta p ~ T play an 
essential role in modifying the dynamics of the p ~ g^T degrees of freedom [^|-|^. Bodeker 
has gone further, showing the role of scatterings between such hard modes by exchange 
of g^T ^ p ^ gT modes There has also been progress in numerically studying the 

effective infrared (IR) theories to determine the size of baryon number violation [p!0| , p!T[] . 

The rate of baryon number violation is characterized by the sphaleron rateQ 



I. INTRODUCTION 
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Our normalization is Di = di + iT^-Af, and Fij = [Di, Dj]. a means a^, and g means g,^ 



Bodeker has shown that in pure Yang-Mills theory, F has the parametric form 



r = .'(log^ + C + 0(l/log)) (^]»'T\ (1.3) 



The constant^C can and recently has been determined analytically |T2|. The leading coeffi- 
cient k' can be obtained numerically by studying a local, UV finite theory, which is precisely 
Langevin dynamics for classical 3 dimensional nonabelian gauge fields: 

Ha = J d'x—F^^Ft^{x) , -g'^ = {D,F,,r = (D x B)^ , (1.5) 

{^ixi, ri)el(x2, rs)) = 2g^T6{xi - X2)6{n - r2)6,,6''' , (1.6) 

Here is the covariant derivative in Langevin time r, F^^j the (nonabelian) magnetic field, 
and ^ Gaussian white noise normalized as shown.Q 

The form of this theory is identical to the nonabelian Ampere's law, but with the current 
replaced with j = crE, as we would expect in a conducting medium; 

Dx B = -aE ( + noise) , (1.7) 

with the Langevin time identified as r = at. Here a is a "color conductivity" which describes 
the current response to very infrared external fields. It can be treated as a constant, to 
(next to) leading order in \og{l/ g), because the mean free path for color changing collisions 
is shorter than the 1/ g'^T scale where F is set, by a factor of \og{l/ g). An explicit expression 
for cr at next to leading log order, found by Arnold and Yaffe, [jl3[, will be presented in Eq. 
(^.23|). is determined in this theory as 



.' = 5^ lim (((^cs(r) - iVcs(0))'-> 

2ti v,t^oo Vt 



This effective theory was studied numerically in [|14|, with the result (see below) that k' ~ 10. 
This effective theory is not sufficient to determine the 0(1/ log) coefficient in Eq. ( |1.3| ), and 
neither will anything we discuss. 



Most of what we have described, and in particular everything in references 0,|^-|T2|,|T^,|T3 
is for the case of pure Yang- Mills theory. This is appropriate if we are interested in temper- 
atures very much higher than the electroweak phase transition (or crossover) temperature, 
because in that case the Higgs boson has a large thermal mass and can be removed pertur- 
batively. However, the main application we are interested in, baryogenesis, requires knowing 
r in a theory with at least one Higgs boson, in the symmetric phase but at or slightly be- 
low the equilibrium temperature Tgq for the electroweak phase transition. It is not clear 



^Actually it is not a constant, but contains a log(log(m£)/(/^T)), see Eq. ( |3.23 ). 

•^Note that the Langevin time r has dimensions of length squared, and recall our scaling convention for A. 
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whether discarding the Higgs physics is justified in this case, and we should rethink both 
the appropriate effective theory, and the numerical determination of F, in this light. 

In this paper we will attempt to fill this gap. First we discuss how important we ex- 
pect Higgs physics to be by considering the thermodynamics of Yang-Mills Higgs fields in 
Section Then we construct an appropriate infrared (IR) effective theory which general- 
izes Bodeker's effective theory in Section |T|. We discuss the numerical implementation in 
Section ^ and present numerical results in |V]. We find that, where the phase transition is 
strong enough to preserve baryon number after its completion, the change in F due to Higgs 
physics is quite a small effect, so the error in using pure Yang-Mills theory is small, of order 
20%. While we work in the minimal standard model (at experimentally excluded values of 
the Higgs mass, to get a strong enough phase transition), we expect the results to hold as 
well in extensions such as the MSSM with a light scalar top quark. This conjecture could 
be tested by simulations in that theory, along the lines of what we do here. 



II. THERMODYNAMIC INFLUENCE OF THE HIGGS 

To get an idea of how important Higgs physics will be for the sphaleron rate, we will try 
to get an idea of how important it is thermodynamically for the infrared transverse gauge 
boson excitations which we expect to be responsible for baryon number violation. As we 
will see, this in fact gives a reasonable estimate for what difference the Higgs will make in 
the sphaleron rate. 

To a very good approximation the thermodynamics of infrared bosonic fields in the hot 
electroweak theory can be described by a three dimensional path integral |TB|. In fact this 



can be understood as a special case of the statement that the IR physics is essentially 
classical, since the three dimensional path integral we arrive at coincides with the partition 
function of the classical bosonic theory. Up to parametrically suppressed corrections^, the 
partition function governing the thermodynamics is |l( 



Z = J VA^VAoV<^exp-H/T , (2.1) 
H = J d'x^^Ft^Ft^ + ^(AAo)^(A,Ao)^ + ^^Ao^S + ^Mo^o 



\a Aa\2 



A. 



+ (A$)^(A$) + m2,(T)$t$ + A(<l>^$)2 + -^AgAS$t<l> . (2.2) 

The couplings A, g"^, \a ~ a^, and \a<s> = fi'^/2 + O(a^) are determined by a matching 
calculation; g"^, A, and the field wave function normalizations correspond to those of the 4-D 
theory at a renormalization point given roughly by T. We will always treat A ~ (7^ for power 
counting purposes, as is appropriate given the renormalization structure of the theory. It is 



''The 3-D theory should be viewed as an IR effective theory for the thermodynamics below the scale T. 
At some level of precision it becomes necessary to include high dimension operators. If we are interested in 
thermodynamics at the length scale then neglecting the high dimension operators causes errors of 

0{{g'^T/T)^) times an additional explicit factor of a because the high dimension operators are radiatively 
induced, leading to an error. Of course, on less infrared scales the effective theory is less accurate. 
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(b) 

FIG. 1. gauge field self-energy insertions generated by the Higgs fields. 

sometimes useful to consider \ <^ and to expand in X/g"^, but if we do so it is implied 
that A is still ^ g^, for instance. Aq is the remnant of the temporal connection and can be 
thought of as an adjoint scalar field. Its mass term is the Debye mass responsible for charge 
screening, and is given by 



mliT) = ^^(4A^c + Nf + 2Ns) - 3D counterterm + 0{g^T'^) , (2.3) 



with Nc = 2 the second Casimir of the group, Nf the number of chiral, fundamental repre- 
sentation fermions, and A^^ the number of fundamental representation complex scalars. In 
the minimal standard model, the 0{g^) piece is = (11/6)5'^T^. 

Because of the Debye mass, the Aq field is heavy and has little influence on the very 
infrared thermodynamics, leading only to a rescaling of the effective gauge coupling by an 
0{g) correction However the same is not generically true for the Higgs field. This is 



because, unlike the Aq field, it has a negative vacuum mass squared which may approximately 
cancel the positive induced thermal mass; 

™?,(r) = m?i + mL, mL<0, m?i ~ g'^^ (2.4) 
(specflcally. = + 8A)r^ i„ a, mSm) . 

At very large temperatures, m|^(T) ~ g^T"^ and its thermodynamic effects are parametrically 
suppressed. However, the electroweak phase transition occurs where the positive thermal 
mass squared and the negative tree one cancel up to 0{g'^T'^) corrections, and it is in this 
regime that we need to know the sphaleron rate. It makes sense, then, to treat m'\j{T) 
parametrically as 0{g*T^), though depending on the strength of the phase transition, it 
may either be "large" or "small" within this parametric range. 

Let us now consider the thermodynamic influence of the Higgs boson on the Yang- 
Mills fields. Since we are only making estimates here, we will be content with one loop 
calculations even where the loop expansion is not very reliable. We can estimate the Higgs 
field's importance by considering Higgs contributions to the gauge field self-energy at soft 
external momentum p ~ g'^T. The sum of diagram (a) in Fig. |I| and any p independent 
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contributions from (b) vanish in any gauge invariance respecting regularization; the scalar 
contribution to the gauge field self-energy is then 



(2.5) 



2 7 (27r)3 V((A;-f)2 + m2,)((A;+f)2 + m2,) (k^ + J 

If we take p <^mH, the result is 

n??(p) = -^ {S,p'-P.p,) , (2.6) 

whereas in the regime where rxiu <^ p, we get 



32b| 

In either case the sign of the effect is such that it reduces the size of infrared gauge field 
excitations; so the presence of a Higgs field should reduce F. 

The calculation is for the symmetric phase, where there is no Higgs condensate. When 
there is a condensate 0o (normalized as 0g = (2$''"$)) then there is an induced mass for 
the gauge bosons of = 5f^0Q/4, which will of course further reduce F. However, we are 
primarily interested right now in the symmetric phase case. 

To decide how important the scalar contribution is, we need to know what momentum 
p is characteristic for the baryon number violating processes we intend to study. Since the 
physics we are after is nonperturbative physics, we expect the answer to be, p such that 
perturbation theory is breaking down. We can estimate what p is necessary by computing 
the contribution of gauge bosons to the self-energy. The result is not gauge invariant, but 
in Landau gauge the analogous contribution from gauge and ghost loops is 

Uijip) from bosons = - piPj^ . (2.8) 

The sign is the opposite, indicating that gauge self-interactions lead to a more rapid onset 
of nonperturbatively large fields. The loop contribution comes on order the tree inverse 
propagator p'^6ij —piPj, and the calculation therefore breaks down completely, for p ~ g^T/3. 
We take this as a fair estimate of the nonperturbative scale, though such an estimate certainly 
cannot be considered accurate to better than about a factor of 2. 

For a heavy Higgs boson, then, we expect a correction of order {g'^T /rriHiT)) / {ASn); but 
even for mniT) ~ we only expect a correction from the Higgs boson to the propagator at 
the relevant momentum to be of order 1/11. This corresponds to a rescaling by ~ (l/H) 
of the length scale where perturbation theory breaks down. However, since the fifth power 
of this length scale enters the sphaleron rate, the importance of the Higgs boson would be 
a little larger, as large as a (10/11)^ ~ 0.6 reduction of F. This analysis also implies that. 



^In Feynman gauge the 11 becomes 14. The answer differs frora that in jl^], for instance, because we have 
left out the (heavy) Aq field, treated as light there. 
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when there is a condensate, it starts to significantly suppress baryon number violation when 
the induced gauge field mass is of order (yf^T/S, which requires 0o — 2(y'T/3. Of course for 
such a small condensate, the perturbative notion of "condensate" is lost and what we are 
doing is unreliable; but this should give us some estimate of how strong the transition needs 
to be before the drop in the sphaleron rate becomes significant. 

In the standard model and in the regime where the phase transition is strong, we can 
estimate the symmetric phase m//(Teq); at leading order in A/(7^, a standard calculation 
from the curvature of the one loop effective potential gives 



Plugging this into Eq. ( |2.6|) gives a self-energy correction of —y A/1 8 (7^, and a correction to 

the sphaleron rate of 



Below we will study the case where \/ g'^ = 0.036, for which the correction is 22%. However, 
this mass is not in the regime mn ^ p for p = g^T/3, so we have computed Eq. ( p.5|) 
numerically for this itlh and p; it yields —0.037. A rough estimate is that we will see F 
reduced to (1 — .037)^ = .83 of the Yang-Mills theory value. This estimate is quite rough, 
but it gives the idea that the symmetric phase rate will be lower, but not much lower, than 
the pure Yang-Mills theory rate. 

To conclude, at the thermodynamic level, one would expect the Higgs field induced 
correction to infrared gauge field behavior, and hence to the sphaleron rate, to be rather 
small in the symmetric phase, and of order 1 for condensates smaller than 0o — (2/3) (7T. 
Naturally, for larger condensates the suppression becomes very substantial, see [0,|18|. 

III. EFFECTIVE THEORY IN THE PRESENCE OF THE HIGGS 

Now we will go about constructing appropriate effective theories for studying baryon 
number violation when there is a light Higgs boson. Our treatment will follow very closely 
that of [^,0,3. Indeed, most of the complications stem from the Yang-Mills sector and have 
been resolved in the literature cited; the Higgs boson will introduce comparatively minor 
new complications. 

We will also leave out the U(l) field in what follows, even though it is light. This is a 
common and probably reasonable approximation. We make it partly because we expect the 
infiuence of the U(l) field to be weak; there is no direct interaction between the SU(2) and 
U(l) gauge fields, and the U(l) physics makes rather small modifications to the thermody- 
namics of the SU(2)-Higgs system |T^. Also, including the U(l) physics would prove to be 
a significant complication, since at the k ~ g^T scale the U(l) fields are overdamped but 
with a k dependent damping which must be treated as a nonlocal effect even in the final 
effective theory. We will not try to address this problem here. 

Since we are not concerned with cases where there is a very large condensate 0o ~ T/g, 
degrees of freedom with momentum k T behave at leading order as massless free fields 




(2.9) 





(2.10) 
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propagating in the background of the IR fields, and the degrees of freedom with k < gT have 
large occupation numbers and can be treated as classical fields. The infrared behavior of 
the theory is described by a classical effective theory in which the ~ T degrees of freedom 
have been analytically integrated out. The resulting equations of motion are similar to those 
of the classical field theory. Defining 

H = Ha + H^ + He + Hu, (3.1) 

(3.2) 

(3.3) 

(3.4) 

(3.5) 

the effective field equations will be the classical equations of motion derived from this Hamil- 
tonian, supplemented by the hard thermal loop effects which arise from the integration over 



Ha = ^ 


1 d'^^^PrM^) ' 


F,,^[D,,D,], 


= 




) + mj^^^^{x) + A(<l>"^$)2(x) 


He = ^ 




Ei = [Dt,Di], 


Hn = ^ 


f d^x^u^u{x) , n 





the heavy modes. The hard thermal loop (HTL) field equations are pO|-p3 

OHa , dH^ rda 



An(x) = -^^^ - mmx) 



9$t(a;) 

= DiDi^x) - {m]j{T) - 2A$^$(a;)) ^{x) , (3.7) 
{D,Wr{x, v) = -v,{D,Wnx, v) + m^ViE^ix) , (3.8) 

(D.E.rix) = (uUn{x) + c.c.) + mj, J ^W'ix, v) . (3.9) 

These equations determine the field evolution, given initial information for A, E, $, 11, and 
W, up to the freedom to choose the time dependent gauge. The last equation is Gauss' Law 
and needs to be applied as a constraint on the initial conditions; it commutes with the other 
equations so it remains valid at later times. 

The only hard thermal loop effect on the soft scalar field is a thermal mass squared 
correction. We discuss this point at some length in Appendix 0, where we demonstrate 
the absence of all other HTL effects involving scalar external lines. We write the correction 
as m^^ (the thermally induced mass squared) and the sum + m^^ = mj^lT). In what 
follows we will absorb this shift in into H<s>. This does not require any modification of 
Eq. ( |3.6| ) because m^j^$^$ is independent of A. 

The HTLs for gauge fields are much more complicated, and are implemented here with 
auxiliary fields, the W fields, which allow them to be written in a local way [^,^. The 
penalty is that the W fields depend on direction v as well as position x. Here f is a unit 
vector, and dfly is an integral over directions (the unit sphere) normalized so / dfl^/Air = 1. 
We normalize W slightly differently than the references, absorbing a factor oi mo into its 
normalization so it enters the Hamiltonian with the same weight as the gauge fields. Its 
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listing in Eq. ( pTT|) would be 



= / d'x— / ^W''{x,v)W''{x,v) . (3.10) 
J 2g^ J ATI- 
It is possible to derive the field equations from this (generalized) Hamiltonian, together with 
Eq. ( p.l|) , but it requires rather nontrivial Lie-Poisson brackets ||2^. Note that the hard 



scalar contributions to the gauge HTL is identical in structure to that from hard gauge and 
fermionic degrees of freedom. This is most easily seen in the kinetic theory derivation of the 



gauge HTL's where one finds that the spin of a hard degree of freedom only matters at 
subleading order in g. 

Eqs. ( p.6D -( ^^ already give an effective theory which is amenable to numerical treat- 



ment along the lines of that given for pure Yang-Mills theory in The added com- 

plication of including the Higgs field is much less than that involved in treating the hard 
thermal loops. However, such a numerical implementation is not ideal, because the system 
of equations presented is not UV finite. In particular the classical gauge and scalar degrees 
of freedom, with dynamics determined by Eqs. ( |3.6| ) and (^.7|), generate UV divergent loop 
corrections which can be considered as extra contributions to the hard thermal loops. In 
a lattice regularization, these extra contributions are finite but grow linearly as the lattice 



spacing is made smaller, and are not rotationally invariant So long as is kept suit- 
ably large this should be a subdominant effect, but it is not clear that it is safe to make mn 
as large as the inverse lattice spacing. Because of these problems it is not easy to extract 
high precision information from such a simulation, and if we expect a fairly small change to 
r due to the Higgs bosons, it may be problematic to isolate it from lattice spacing effects. 
Rather, we will follow Bodeker and integrate out more degrees of freedom to construct 
what ultimately proves to be a simpler and cleaner effective theory (though no longer valid 
to corrections parametrically suppressed by a full power of g), which will prove a better 
test-bed for studying the importance of Higgs physics. 

With this in mind we consider integrating over the gT scale, down to some intermediate 
scale /i -C gT. We are also removing physics with a frequency scale u ~ gT; however we 
cannot directly consider only degrees of freedom with <^ A; ^ gT; for now u is permitted 
to be as large as fi. 

Here we will only treat the case in which there is not a Higgs boson condensate of order 
00 ~ f^/g- This either restricts how deeply into the broken phase our analysis remains valid, 
or how small we are permitted to make the scale fi which we integrate down to. We do 
not consider this restriction problematic because we are mostly interested in the symmetric 
phase, and when there is a large condensate we will need to use other tools to determine the 



rate anyway. We have already considered the broken phase problem in ^8 . 

Under this assumption, the behavior of gauge bosons with k ^ ^ is not significantly 
changed by a Higgs condensate. It is also not significantly changed by the background of 
IR Higgs fluctuations, as can be verified by a loopwise analysis like the one in the last 
section; and the effects of hard Higgs excitations are already included by the HTL effective 
action. Therefore the integration over these degrees of freedom proceeds as it does in the 
case without Higgs fields. The integration over gauge and W field degrees of freedom is 
very nontrivial and has been treated at length by Bodeker |^J^ and by Arnold, Son, and 
Yaffe 0JT^ . They show that a collision term is induced for the W fields, together with noise 
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required by the fluctuation dissipation theorem. The integration over the gT scale Higgs 
fields is much simpler. The Higgs field equation of motion, Eq. ( p.7|) , is the same below the 
gT scale as it is up to the T scale. Provided that we keep /i ^ mH(T), the Higgs field still 
undergoes free relativistic propagation in the soft gauge and Higgs field background up to 
corrections 0{g'^T/fi). Such propagation is precisely what generated the hard thermal loops. 
Standard power counting shows that contributions from the Higgs fields with momentum 
k > fi to more IR degrees of freedom are suppressed by g^T/fi except for UV divergent 



contributions. But, as shown in [Q, the structure of UV divergences in the classical field 
theory coincides exactly with the hard thermal loops. Hence the k > fi Higgs degrees of 
freedom induce HTL's and parametrically suppressed additional effects. The HTL effects 
from the gT Higgs fields are smaller by a factor of g than those from T physics, so we can 
neglect them next to the HTL effects already being included. 

It is worth remarking why physically no collision integral is induced when we integrate 
out scalars with momentum k ~ gT. The total scattering rate for hard modes by exchange 
of a soft gauge particle is ~ g'^T \og{l / g) , with the logarithm arising from the momentum 
region gT to g^T. This large collision rate originates from the IR singular s^/t^ matrix 
elements in gauge boson exchange. When we integrate out gT bosons we must explicitly 
include their contribution to scatterings via a collision integral. On the other hand, no 
scattering process between hard modes which is mediated by single scalar exchange has a 
scattering rate greater than ~ g^TloglT/rriH)^ Such rare scatterings are not important at 
leading order because we are eventually interested in physics at the 1/g'^T length scale. 

On length scales more infrared than 1/gT, Eqs. (|3.6| ) and ( p.8|) describe overdamped 
evolution for the gauge fields. Further, the W field equation of motion is linear in its slowly 
varying source oc Wj-Ej. Therefore it is permissible at leading order in fi/gT to drop both the 
DtE term in Eq. ( |3.6D and the DfW term in Eq. leading to an effective theory 

for the /i ^ gT scale physics, [] 

m^J -V.-W {x,v) = -g , (3.11) 

AH(x)^-^, (3.12) 

v,{D,WY{x, v) = m^v,E^{x) - J ^C,,W{x, v') + C(x, v) , (3.13) 
TTiD f^Wix^v) = 0. (3.14) 



The new features are the collision integral C„„/ , whose form and value is discussed in [|^^,|T2[ , 
and the noise (. At leading order in \og{gT/fi) the collision integral is 

C^v' ^ 1 \ Ss^iv - V ) , ^-—A log • 3.15 



^This can be verified by looking at the matrix elements of tree level 2 — > 2 scattering processes mediated 
by a scalar, none of which contain s^/i^ type terms. At worst they go as s^/tu and are log IR divergent. 

'^For a discussion of the last equation, see Pq]. 
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The noise is Gaussian and white with two-point correlator 



V, t)C\y, v', t')) = 2g'TC,,,6'''6{x - y)5{t - t') . (3.16) 

Because of our nonstandard W field normalization, the normalization of the noise differs 
from the references. Note that the size of the collision integral is parametrically Cy^i ~ 
g^T\og{gT/^^). 

At the k ~ g^T scale this effective theory has two natural time scales. First there is the 
time scale ~ 1/g'^T, on which Eq. ( |3.12| ) allows the Higgs fields to evolve. There is also a 
scale, ~ 1 / g'^T log{l / g) , on which the overdamped gauge and W fields evolve. This reflects 
the fact that the gauge HTL's include Landau damping and lead to overdamped evolution, 
while the Higgs HTL's are just a mass correction and do not induce any g'^T'^ size damping, 
so the Higgs field is not overdamped. 

If we view the system on length scales ~ l/g^T and on the Higgs time scale ~ 1/g'^T, 
the gauge fields look "frozen" up to parametrically suppressed effects. On these time scales 
the Higgs field evolution is determined by Eq. ( |3.12| ) with the gauge field background 
frozen. What is the behavior of such a system? The Higgs propagates on an inhomogeneous 
background connection, interacting with itself via the nonlinear quartic coupling term. There 
is a widely (though not universally) believed conjecture that the evolution of a scalar field 
theory with quartic self-coupling in 3 dimensions should be ergodic, in which case it will 
randomize itself on the 1/g'^T time scale (since there is no other available time scale for 
its evolution). This view is supported, for instance, by the calculation of the damping rate 
of a scalar at rest, in a fiat connection but with a quartic interaction; the damping rate 
is parametrically ~ X^T'^/rriH ~ g^T, and the damping arises primarily from degrees of 
freedom with k ~ rriH We expect that the inhomogeneous connection should only 

make the randomization of the scalar field more efficient. In particular we speculate that 
the spectrum of the operator for a typical 3-D gauge field background exhibits Anderson 
localization at all frequencies. 

On the 1/g'^T time scale on which the gauge field evolves, the Higgs field will thoroughly 
explore its fixed connection thermal ensemble. Therefore, on the time scale on which A 
evolves, it will see a thermodynamic average of the possible Higgs field configurations. We 
emphasize that we are relying here on the conjectured ergodicity of 3-D scalar 0^ theory. It is 
also not clear that our treatment will remain true very near the endpoint of the electroweak 
phase transition, where the Higgs field correlation length grows to be ^ 1/g'^T, because 
in this regime it may take much longer for the Higgs field to explore its fixed gauge field 
ensemble. We will exclude that regime from consideration, although it is not clear to us 
that the effective theory we will derive cannot be used there as well. 

When our assumption is valid, we should average the $ dependent part of Eq. ( ^.11[ ) 
over the Higgs thermal ensemble. 

The RHS of Eq. ( p.ll| ) can be understood as g^T times a variation of a 3-D effective action 
describing the thermodynamics of the gauge fields. The sole modification from the inclusion 
of Higgs fields is that the effective action should include, besides the gauge part, a nonlocal 
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piece arising from integrating over the Higgs fields. Hence tlie effective theory describing 
infrared gauge bosons is 



mo 



Air 



v,W\x,v) = -g'T 



2^ d H^siA) 



dAf(x) T 



— log j P$ exp 



Ha + 



T 



Vi{DiWY{x, v) = mj,v,E1{E) - j ^C^^W^x, v') + v 



(3.18) 

(3.19) 
(3.20) 



This is identical with the "theory 2" of reference ||13| . |1^ , except for the Higgs field additions 

to i/eff- 

This theory is probably not well suited to numerical study. However, as Bodeker has 
shown, it has fairly simple behavior when studied at the length scale 0. In this 

regime the collision integral dominates over the derivative term for the W fieldQ, since 
C ~ g^T \og{gT / fi) with fi ~ g^T, while the derivative term is f ■ -D ~ g^T; the collision 
term is therefore bigger by ~ log{l/g). Roughly speaking, one may solve for W in terms of E 
in Eq. ( |3.2UD and plug it into Eq. ( p. 181 ), yielding a local expression of form aE = —dH/dA. 

At leading log order the Yang-Mills theory argument goes over directly to the case 
including a Higgs field, since it depends on manipulations of Eq. (|3.20| ) only and this does 
not include the Higgs field. The next to leading log calculation of ||12[ also still holds, 
provided there is no Higgs condensate of size 0o ~ gT/\og(l/g). We discuss this point 
at more length in an appendix, where we show how the addition of the Higgs field does 
not modify the calculation presented in |jl2i. Intuitively, the reason the Higgs field does 



not affect a at next to leading log order is as follows. The conductivity depends on the 
efficiency of collisions. The collisions are mediated by gauge excitations with momenta 
gT ^ k > g'^T \og{l / g) . At the low end of this range the Higgs fields modify the gauge 
field thermodynamics by 0(1/ log(l/5')), see Section and this part of the range gives a 
contribution down by l/\og{l/g) compared to the complete collision integral. Hence the 
influence of the Higgs field is suppressed by two powers of log. 
Finally, we arrive at the effective theory 



aEf(x) 



d H,s{A) 



dAf{x) T 
t)ej{y, t')) = 2g^Ta5,,5a,5\x - y)5{t - t' 
3 N,g^T 



-1 



—7, 



7 



An 



In ^ + 3.041 

7 



(3.21) 
(3.22) 
(3.23) 



Here i/efr is as in Eq. (|3.19| ) and the expression for is from [|T3| with a particularly nice 
choice for the renormalization scale /i. The sole modification the inclusion of the Higgs fields 
has made is in the form of -ffefr- Again we emphasize that the the derivation has assumed 
that the Higgs is light and that there is not a large Higgs condensate. If there is a large 



^What follows is a gross oversimplification of the argument, see [^1). 



11 



Higgs condensate 0o ~ gTl\og{l/g) then the value of a changes, and if 0o ~ ^ then the 
form of the effective theory changes as well. We will not discuss the latter case here. Also 
note that our effective theory is only valid for studying gauge field correlators; it cannot tell 
us much about unequal time Higgs field correlators. 



IV. NUMERICS 



Bodeker's effective theory is Eq. ( |3.21| ), but with i^efr replaced with Ha- It makes an 



excellent starting point for numerical investigation of F in pure Yang-Mills theory for two 
reasons: 

1. it is local, and 

2. it is UV finite. 

Neither is true for the effective theory we have derived. This potentially makes its study 
much more problematic than Bodeker's effective theory. 

In practice the second problem is not a substantial one. The theory is not UV finite 
because H^s involves the 3 dimensional path integral for a scalar field, which contains linear 
and logarithmic mass divergences arising from one and two loop graphs. However, the 
path integral is still super-renormalizable, and the UV infinities are purely local and can 



be absorbed with a mass counterterm. Their value is known [^^, and as we will discuss 
below, we could actually proceed even if they were not. 

The real problem with implementing the effective theory is its nonlocality, which comes 
about because of the path integral in the expression for i^efr, see Eq. ( p.l9| ). The solution 
is to think about how we would carry out such a path integral numerically. To evaluate Eq. 
( p.l9| ) numerically, we would perform the path integral by Monte-Carlo, for instance by a 
Langevin equation, 

{f(i,r^)f'(i',T;)> = 2T lS{x - x')S{Tt - r;) , (4.2) 

Here the 1 reminds us to make the noise diagonal in the components of the Higgs field. 
This Langevin equation must be evolved, and the result averaged, at every time step in the 
Langevin dynamics of the gauge fields. 

This suggests that Eq. ( p.21| ) can be replaced with the rj ^ (yo limit of the following 
system of equations: 

aEtix) = -g'^^^(HAiA) + H^iA, $)) + (x, t) , (4.3) 
d 

aDMx) = -r]-^^^H^{A, $) + ^^{x, t) , 
t)^Ux', t')) = 2r]aT 16{x - x')6{t - t') . 

This is Langevin evolution but with a different rate for the evolution of the two fields, A 
and $. In the limit that the Higgs field evolution is made infinitely fast, which is 1] oo, 
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the Higgs field evolution will perform the path integral in a much shorter time scale than 
the gauge fields evolve in, and we recover the desired equations of motion, Eq. (p.21|) . 



Eq. ( |4.3| ) at a finite value of r] does make a good starting point for numerical work. 
The large rj limit must then be taken numerically. It is also convenient to rescale time to a 
Langevin time, r = at. Note that the dimensions of r are those of a length squared. 

Our nonperturbative regularization for Eq. ( [4.3|) will be the lattice. Since all of the 
numerical tools we use exist in the literature, we will only present the relevant references 
here rather than give complete details. The lattice discretization is standard, see |3^. Our 
topological lattice definition of iVcs is the same as in |^]. The relation between the couplings 



of the lattice and continuum systems has been worked out in and we use the expressions 
there. These relations match all thermodynamic quantities at 0(a), leaving O(a^) errors, 
with two exceptions. We do not know the full 0{a) match for an additive piece of the 0^ 
operator insertion or for m'jj{T). If we were interested in determining F at a particular, 
fixed value of m'\j{T) then this could pose a problem. However, what we want is F when 
m'\j{T) is a fixed distance from the value which gives phase equilibrium, m1^. To determine 
this we do not need the absolute normalization of m'\j{T) to 0(a), but we will need to find 
TTigq numerically. 

Now we briefly discuss algorithm. To determine m^q we use the same algorithm and 



multicanonical techniques as [g^]. To perform Eq. (|0|), we may use whatever Higgs update 
we choose, provided that it is stochastic and that we take the rj —>■ oo limit. We choose 
a mixture of heat bath Higgs updates and the x-y over-relaxation algorithm of [^. For 



the gauge field update we should use Langevin dynamics or any other strictly dissipative 
dynamics. We choose heat bath. For either algorithm, the effective infrared dynamics should 
be of the Langevin form, and the difference between the algorithms will be a radiative 
rescaling of the (Langevin) time scale and high dimension corrections, which first appear at 
O(a^)0. The radiative rescaling of the time scale is discussed at length in Appendix A of 
. That paper presents an analytic calculation of the lattice to continuum match for the 
Langevin algorithm, and the Langevin to heat-bath rescaling is found numerically. Within 
error it is equal to the coefficient which appears in the 0(a) thermodynamical match 
detailed in . We will assume that this numerically determined relation holds analjd^ically. 

V. RESULTS 

With a numerical implementation of the effective theory in hand we will address two 
questions. First, how large is F in the symmetric phase, at ml^, for parameters where the 
phase transition is strong enough to preserve baryon number after its completion? Second, 
how does F vary as we go through the analytic crossover present for large vacuum Higgs 
masses? 

To address the first question we consider Yang-Mills Higgs theory at X/g'^ = 0.036, which 
is roughly the value at which the phase transition is barely strong enough to preserve baryon 



^The question of what modifications may occur in the effective IR behavior has been addressed at some 
length by Arnold and Yaffe, see particularly Appendix A of |l5| 
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FIG. 2. Log of probability distribution for </)^ = 2$^$ (MS renormalization point JI = g'^T) for cubic 
boxes of size (left) and 32'^ (right), at lattice spacing A/dg^T {(3 = 9). Inverting the y axis roughly gives 
the free energy dependence on 0^. The left peak is the symmetric phase and the right peak is the broken 
phase; the dip between them is a free energy barrier separating the phases, which grows as the volume is 
increased (it is not free energy per volume). The asymmetric shape is typical for this order parameter. 



number cosmologically after its completion.[^ We determine at three lattice spacings, 
a = 4/7g'^T, a = 4/9g'^T, and a = 4/12g'^T (which in the notation which is sometimes 
customary [Rl] are j3 = 7, 9, and 12), by multicanonical Monte-Carlo. The probability 
distributions for 0^ at for two lattice volumes and a = 4/9g'^T are shown in Fig. ^. 
First we consider the t] ^ oo limit at fixed lattice spacing a = 4/9g'^T in a cubic volume 
with periodic boundary conditions, for L = 14.2/ g'^T (32 sites on a side). This volume 
is abundantly large enough to see the large volume value of F ||3l| , and is large enough that 
strong metastability prevents tunneling to the broken phase at ml^, where we work. Table 
I presents k', defined in Eq. (|1.3| ), for various values of rj, lattice spacing, and mjj. The rj 
dependence is weak and statistically compatible with zero. 



a X g'^T 


V 




volume X Langevin time 


k' 


4/9 


2.5 





32^ X 56400 


8.45 ± .24 


4/9 


5 





32^ X 56400 


8.20 ± .24 


4/9 


10 





32^ X 56400 


8.08 ± .23 


4/9 


5 


.0093 


32^ X 34300 


8.63 ± .32 


4/9 


5 


-.0093 


32^ X 30900 


8.50 ± .33 


4/7 


10 





24^ X 57400 


8.13 ± .37 


4/12 


5 





40^ X 98600 


7.83 ± .27 


4/9 


Yang 


Mills 


243 X 489400 


9.90 ± .13 


4/10 


Yang 


Mills 


40^ X 84900 


10.00 ± .23 



The value we use is slightly smaller than that found in 113] because here we do not include the U(l) field. 
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5000 10' 



Langevin time (lattice units) 
FIG. 3. Left: probability distribution on supercooling, mf^ = — .0093(7'*T^. Right: Langevin time 
history of A^cs starting in the symmetric phase, under such supercooling. After a tunneling event to the 
broken phase, A^cs stopped diffusing. 



TABLE L Chern-Simons number diffusion F in the symmetric phase near the equilibrium temperature 
with A/g^ — .036, for different lattice spacings, values of rj, and Higgs masses. The 4- volume is expressed in 
lattice units; to convert to physical units multiply lengths by a and Langevin time by aa? . Except for the 
pure Yang-Mills data, all values of k' are statistically compatible; no lattice or -q dependence is statistically 
significant. 

Table | also shows that, as expected, there is almost no lattice spacing dependence in 
r. The 0(a) match for the thermodynamic quantities and time scales is essential here, 
see |T^. Finally, since in the usual electroweak baryogenesis scenario the symmetric phase 
undergoes supercooling, we should study V in the supercooled symmetric phase. Table 
I shows that the maximum supercooling compatible with strong metastability is still not 
enough to significantly change F. We cannot increase the supercooling beyond what was 
used because the lattice system will nucleate to the broken phase. Indeed, one run used for 
the table ended with a nucleation to the broken phase, as is seen clearly from the time history 
for A^cs in that run, shown in Fig. |^. Naturally, the broken phase part of the evolution was 
not used in the analysis. Combining all the figures in the table, since all are statistically 
compatible and no trend (lattice spacing, supercooling, or rj) is statistically significant, we 
get 



F ~ 

symm — 



8.24 ±0.10 




log^ + da^T\ (5.1) 



for X/g"^ = 0.036. We have computed the pure Yang-Mills theory value at the same lattice 
spacing with comparable precision to facilitate comparison; the symmetric phase value is 
0.832 ± .015 of the Yang-Mills theory value. (The Yang-Mills theory value found here is 
smaller than that quoted in |jl^, where we found k' = 10.8 ± 0.7. This is a statistical 



fiuctuation in the data there. The table includes a higher statistics redetermination of k' at 
the same parameters used in that paper, with a result in agreement with the [3 = 9 value 
found here, and statistically compatible at Icr with the result determined in [H.) 
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m|(T)/g^T2 rnf,(T)/g*T2 
FIG. 4. Left: (j)^ susceptibility; Right: (f^ susceptibility (peak with small errors) and sphaleron rate 
(larger errors, square plotting symbols) scaled to their maxima, when there is a smooth crossover. m^{T) 
plotted is the 3D theory value, MS renormalized with JI = g^T. 



Of course the above is academic because X/g'^ = 0.036 converts to a vacuum Higgs 
mass rriH = 44GeV at tree level. Including the large radiative top Yukawa corrections, 
this X/g"^ does not correspond to any physical Higgs mass |jl6[. In the standard model, all 



experimentally allowed Higgs masses fail to provide an electroweak phase transition at all. 
Rather, there is an analytic crossover [R4|. Although baryogenesis is probably impossible in 



that setting, it would still be interesting to see how T varies as we go through the crossover. 
Does it turn on suddenly or gradually, and does its turn-on point coincide with the peak 
in the susceptibility for 0^? We answer this question in Fig. ^, which shows how the (j)'^ 
susceptibility 



V 




(5.2) 



and the sphaleron rate F vary with m'\j{T). The data are for X/ g"^ = 5/16, corresponding 
at tree level to a physical vacuum mn = 130GeV, and were taken for a 32'^ box with lattice 
spacing a = l/2g^T. The sphaleron rate data were all taken using rj = 20. We see that the 
switch-on of F, though smooth, is fairly rapid and occurs at slightly lower m|^(T) than the 
peak susceptibility; that is, it is when conditions are a little more "broken phase-like" than 
when the susceptibility peaks. The sphaleron rate F proves a rather good order parameter 
to distinguish Higgs-like and symmetric-like phases. Although, like any order parameter 
must, it shows smooth behavior, the range where it is far from both its broken phase value 
k' ~ and its symmetric phase value k' > 5 is quite narrow, roughly as narrow as the peak 
in the 0^ susceptibility. However, we do not view the sphaleron rate as competitive with the 
susceptibility as a probe of where the crossover occurs. The main reason is that it is possible 
to perform a multicanonical reweighting which allows the susceptibility to be scanned in a 
wide range of m'jj (the data in the figure come from a single numerical run), and the statistics 
for the susceptibility improve more quickly. For instance, the Nqs diffusion data in Fig. 
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took more CPU time than the susceptibihty data, but are substantially "dirtier". (The 
data sets are plotted with their error bars. The A'^cs diffusion errors are each independent, 
but there is very large cross-correlation in errors of neighboring points for the susceptibility 
because they were all computed from one data set.) 

VI. CONCLUSION 

Assuming that classical 3-D scalar theory is ergodic, the addition of a light Higgs de- 
gree of freedom replaces Bodeker's effective theory for the evolution of infrared gauge fields 
with a slightly more complicated equation, Eq. ( |3.21| ). The sole change is in the thermo- 
dynamic potential for the gauge fields. Including the Higgs field makes this thermodynamic 
potential nonlocal. However, the effective theory is still a useful starting point for numerical 
work because we can use the limit of a sequence of local effective theories, namely, Langevin 
(or heat bath) evolution for the gauge fields and for the Higgs fields, but with much faster 
Langevin evolution for the Higgs degrees of freedom. 

When the electroweak phase transition is strong, the sphaleron rate in the symmetric 
phase is reduced by around 20% from its Yang-Mills theory value, roughly in accord with 
an estimate based on the thermodynamics. (The thermodynamic estimate for small \/ g"^ is 

(5r/r) ~ 1.2y^A/(7^, see Eq. (|2.10|) .) When there is no electroweak phase transition, but 
an analytic crossover, the sphaleron rate changes rather quickly from its symmetric phase 
value to nearly zero, roughly at the same value of m^jiT) where the 0^ susceptibility peaks. 
The crossover region is of about the same width as the peak in the 0^ susceptibility and is 
displaced to slightly lower rn?^{T). 

In this paper we have only studied the standard model, either for parameters which are 
ruled out experimentally or for which the electroweak phase transition cannot provide for 
baryogenesis. However it is fairly simple to see how to extend the work to more viable 
models like the MSSM (minimal supersymmetric standard model). In that case, the SU(2) 
and SU(3) gauge fields would each evolve under Langevin dynamics, with a thermodynamic 
potential arising from integrating over the Higgs and scalar top fields. However, from our 
results with a strong phase transition, it should be clear that, in all cases where the phase 
transition is strongly first order, F in the symmetric phase will be very close to its Yang-Mills 
theory value. This is because we have found that the suppression of F corresponds well to 
what we expect thermodynamically; and when the electroweak phase transition is strong 
in the MSSM, the symmetric phase Higgs mass is larger than in the standard model case. 
Therefore the SU(2) thermodynamics in the symmetric phase is closer to Yang-Mills theory 
in the MSSM than in the standard model, see Eq. ( |2.6| ). Hence it is almost certainly true 
that, in the MSSM and when the phase transition is strong, the symmetric phase sphaleron 
rate is lower than but within 20% of the pure Yang-Mills theory value. In practice this 
means we can continue to quote the Yang-Mills theory result for all symmetric phase cases 
of physical interest, with modest error. 
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APPENDIX A: HIGGS FIELDS AND HARD THERMAL LOOPS 

In this appendix we show that the only hard thermal loop needed for the bosonic effective 
theory, which has soft external scalar field lines, is the scalar mass correction. 

We need only consider HTL's with bosonic external lines, since the appropriate IR ef- 
fective theory we seek is bosonic. Hence we need in general to consider all diagrams with 
{us > 0,even) scalar and Ug gauge boson external lines. We will use repeatedly the power 
counting rules derived in [Q, which we repeat here for the reader's convenience. A hard 
thermal loop is always an O(T^) contribution from the ~ T momentum region of a one 
loop integral with soft {gT) external lines. To determine the largest power of T possible 
from a diagram, follow these rules: 

1. the loop integration contributes T^; 

2. the first propagator times the Matsubara sum contributes 1/T; 

3. every additional propagator contributes 1/(PT); 

4. powers of in the numerator (from 3 point gauge vertices or fermionic propagators) 
each contribute T; 

5. when there are two or more propagators, and all propagators are either bosonic or 
fermionic, there is an extra P/T suppression. 

If the result has a weaker power of T than the diagram does not contribute an HTL; if 
it is the diagram will unless there is some cancellation. 

Since we are only concerned with HTL's with bosonic external legs, then either all the 
propagators in the loop will be bosons or all will be fermions. Any diagram with more 
than 1 propagator will get a P/T from rule 5. and this will exclude diagrams with 4 point 
vertices (4 gluon, 2 gluon and 2 scalar, or 4 scalar) except for tadpoles, just as is the case 
in Yang-Mills theory. This leaves all diagrams of the general form of diagrams (c) and (d) 
in Figure ^. 

The scalar self-energy diagrams are considered in |^ (fermion loop) and (boson 



loops). The result is that the self-energy is a momentum independent mass correction. In 
particular it is useful to reproduce the argument for the gauge HTL from diagram (a) of 
Figure |. In Feynman gauge this diagram contributes (following the notation of that paper, 
where capital letters are 4-vectors and lower case are spatial or temporal (0 subscript) 
components) 
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(b) 

(c) (d) 

FIG. 5. Diagrams considered in Appendix Solid lines are scalars, wiggly lines are gauge bosons, 
double lines are fermions. Diagrams (a) and (b) are a self-energy and a vertex correction considered in the 
text. Diagrams (c) and (d) are generic diagrams with external gauge and scalar lines and either a boson 
(alternating Higgs and gauge, diagram (c)) or a fermion (diagram (d)) in the loop. As discussed, (b), and 
all diagrams of type (c) and (d), do not give rise to hard thermal loops. 



2N 



(fk {2P + KY{2P + KYg^, 



2N 



ko 



(27r)3 



(27r) 



2{K + Pf - + 2P^ 
K^{K + Py 



(Al) 



For the first term in the square brackets, the {K + P^ in the numerator cancels a propagator 
in the denominator, leaving one propagator and no powers of K in the numerator. This gives 
an HTL effect only because there is a single propagator left, so 5. above is not invoked. The 
result is a momentum independent tadpole. The second term is similar. The last term has 
two propagators and no power of K in the numerator; it is 0(T°) and does not give an HTL. 

Similarly, diagram (b) gives, in Feynman gauge and labeling the incoming scalar momen- 
tum Pi, the incoming gauge momentum P2, and the outgoing scalar momentum P3 = P1+P2, 



oc 



ko 



d^k {K + 2Pi) ■ {K + 2P^){2K + Pi + P^Y 



{2^y 



K^K 



PiYiK + PsY 



(A2) 



It is convenient to rewrite 



{K + 2Pi) ■ {K + 2P3) = {K + Pif + {K + P^Y -K^ + (4Pi • P3 - P^ - P^) . 

(A3) 

Each K"^ type term will cancel a propagator, leaving two propagators and one K'^ in the 
numerator. By the power counting rules this contributes at most T^{1/T){1/PT)T{P/T) oc 

and does not give a hard thermal loop. The P^ term is 0(T°), so it does not either. 

This example illustrates why diagrams of form (c) and (d), with [ug > 0,even) scalars 
and Ug gauge bosons, (n<j + Ug) > 2, do not give hard thermal loops. Such a diagram 
has {us + rig) > 2 propagators. It also has + Ug factors of in the numerator, either 
all from gauge vertices (bosonic loops) or all from fermionic propagators (fermionic loops). 
However, there are only external Lorentz indices. In Feynman gauge, the diagram will 
be a sum of terms, with at least Ug of the contracted against each other in each term. 
Hence, each term from each diagram contains in its numerator, in Feynman gauge, a term 
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of form {K + Pi) ■ {K + Pj), with Pi, Pj some linear combinations of external momenta and 
K the loop momentum. This can always be written as (sign) (if + PiY + 0{P'^), with 
each {K + P^) a momentum on some propagator. Consider one of the {K + PiY terms. 
Cancelling the {K + PiY against the appropriate propagator leaves (n^ + — 2) factors 
in the numerator but {ug + ng — l) > 1 propagators. The power counting then gives at most 
T3(l/T)(l/PT)"=+"f-2T"»+"9-2(p/T) oc T\ too small to contribute an HTL, even without 
further cancellations. Meanwhile, the O(P^) term has [ns + Ug) propagators but {ns + ng — 2) 
factors; it is at most p2T3(l/T)(l/PT)"^+"«-^T"^+"«-2(P/T) oc T° and also gives no 
HTL. Though we worked in Feynman gauge, we expect the hard thermal loops to be gauge 
invariant; if they vanish except a momentum independent scalar mass in Feynman gauge, 
they should in any gauge. 

Hence the only bosonic HTL with scalar external lines is the gauge field mass insertion. 
Note that similar reasoning to what we give above also rules out HTL's with both scalar 
and fermionic external lines; but this is unimportant for the purposes of this paper. 



APPENDIX B: HIGGS FIELDS AND NEXT TO LEADING LOG 

CONDUCTIVITY 



In this appendix we briefly explain how the calculation presented in |]T2[ is modified by 
the presence of Higgs fields, and why the modification does not change the conductivity of 
the final effective theory at next to leading log order. The purpose of the appendix is to 
outline the argument, at times we will be sloppy with notations and with terms which are 
not relevant and have been discussed at much more length in [13,13,1^1. 

The relevant part of the calculation in ||12[ is the match between their "theory 2" and 
"theory 3", only now "theory 2" is given by Eqs. ( |3.18[ - pl20D , which for convenience we 
repeat here: 



T 



dAf(x) T 



— log J exp 



Ha + Hcf 



Vi{DiWr{x,v) = mnViE^{E) 



T 



(Bl) 
(B2) 
(B3) 



It is possible, at least formally, to invert Eq. 
( [Bl| ) then becomes 



to solve for W in terms of E. Eq. 



(B4) 



where cr{D) is the nonlocal operator resulting from inverting Eq. (|B3| ); it is discussed in |]T2 
(' is a noise with correlator {('(') = 2Ta{D). (Starting here we will aggressively suppress 
indices where we feel the meaning is clear.) 

This Langevin equation only has a path integral representation if we are willing to accept 
{6HeQ/6Ay inside the action of a path integral. It is not clear how to form a perturbation 
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theory for such a path integral. However, if we rewrite Hes as we did in Section then it 
will become possible. The effective theory of interest is the rj oo limit of 



'OH A dH. 



t)f (x', t')) = 2r]T16{x - x')6{t - t') 



(B5) 

(B6) 
(B7) 



To write a path integral expression for this, we follow the standard trick of writing the 
Langevin equation and the average over the noise distribution as a path integral. 



X 



6 (C' + a{D)E + g^T [6{Ha + H^)/6A]) 5 + A$ + [6H^/6^ 



(B8) 



times Jacobians for each field, which are not important, see [T^,[T2|. Here the path integral 
performs the average over the noise, and the delta functions enforce the Langevin equations. 
Enforcing the delta function does the integrals over each noise, giving a path integral 



j VAV^exp{-L) 



L = {a{D)E + g^T 



5{Ha + H^) 



5A 



1 



5JU 



AgMD)T 
6H^ 



a{D)E + g^T 



5{Ha + H^) 
6A 



(B9) 



Here we have not written the contributions from the Jacobians or from an extra regulation 
dependent term arising because of the nonlocality of cr, which ref. calls These 
terms, and the reasons they can be dropped, are discussed at some length in [T5|JT^ . 

Arnold and Yaffe have shown that a in Eq. ( |3.21D can be determined at next to leading 
log order by computing the = 0, A; — limit of the one loop. Coulomb gauge AqAq 
self-energy for the theory with cr(D), and finding what constant value of a is needed to get 
the same one loop result in the final effective theory. To account for Higgs contributions in 
this calculation, we must find all one loop Higgs field self-energy corrections to the Aq field, 
and must determine whether their contributions survive in the 1] oo limit. 

First we have to find what vertices couple $ to Aq. The only terms in the Lagrangian 
which contain Aq at all are those which contain E or Df. They are 



(BIO) 



plus terms odd in Df, which are of form 



2 ' SA. 




+ CC 



(Bll) 
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The first term here is proportional to OiHa- The second is proportional to dtHi^. Both are 
total spacetime derivatives, which can be converted to boundary terms at asymptotically 
early and late times and therefore neglected. Hence we only need to consider vertices which 
arise from the terms in Eq. (|B1CI| ). These allow two new diagrams not considered in [jl2 



which are the same as (a) and (b) in Figure |I| if the wavy lines are now taken to be Aq 
propagators. 

To show that these diagrams vanish in the small rj limit, we need only count powers of 
1] in the vertices and propagators. The Ao$^ vertex carries a factor of lj/t] and the Aq^^ 
vertex carries a factor of l/i]. The $ propagator is 

($$) = — -i — — . (B12) 
The contribution from diagram (a) then scales as 

(a) oc - / dujSk ~ r/-^/^ , (B13) 

?7 J 10^ jri + rik^ 

while diagram (b) behaves as 

and we see that both vanish when we take rj ^ oo. Therefore there are no new contributions 
at one loop from Higgs bosons, and the next to leading order calculation of a is unaffected. 
This will not be true in the case were there is a Higgs condensate of magnitude (po ~ 
gT/log{l/g), because in that case the external Higgs field insertions on the A field lines in 



the one loop diagrams considered in [|T2| will change the gauge field propagators by order 1. 
The Higgs field will also be important beyond one loop, in the construction of an effective 
theory which can determine the 0(1/ log) term in Eq. ( |1.3D ; but we will not investigate that 
problem here. 
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